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ABSTRACT 

The  optical  properties  of  vacuum-deposited  thin  films  were  studied 
in  situ  during  deposition  and  in  air  after  deposition  at  the  Ecole  Nationak 
Superieure  de  Physique,  Marseille,  and  at  the  Optical  Sciences  Center, 
University  of  Arizona,  with  the  objective  of  a  better  understanding  of 
aspects  of  their  behaviour  determined  largely  by  their  microstructure. 

The  basic  measurements  that  were  made  in  both  institutes  were  similar,  T 
and  R  after  deposition  and  T  during  deposition,  both  as  functions  of 
wavelength  in  the  visible  and  near  infrared.  The  techniques  for  reducing 
the  measurements  to  values  of  n,  k  and  of  inhomogeneity  differed  in 
certain  respects  but  the  results  were  similar. 
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INTRODUCTION 

The  purposes  of  this  research  were  twofold.  First  it  was  intended  to 
study  vacuum-deposited  thin  films  by  measurement  of  their  optical 
properties  especially  those  that  would  throw  light  on  their 
microstructure.  Second  it  was  hoped  that  the  research  could  be  used  as  a 
vehicle  to  encourage  collaboration  between  the  group  In  Marseille  and  a 
similar  research  group  at  the  Optical  Sciences  Center  in  Arizona.  This 
report  describes  the  progress  that  has  been  made  In  both  these  areas 
during  the  year  of  operation. 

.The  scientific  objectives  were: 

(a)  The  measurement  of  the  optical  constants  of  thin-films  in  air 
after  deposition. 

Cb)  The  measurement  of  the  optical  constants  of  thin  films  in  situ. 

fc)  The  study  of  the  differences  in  the  measurements  under  (a)  and 
Cb). 

(d)  The  study  of  aspects  of  the  thin-film  deposition  process 

affecting  the  film  performance  as  measured  under  fa)  and^b). 

te)  The  measurement  of  optical  scattering  from  thin  films  and  its 

relation  with  film  structure.., 

/ 

In  one  year  of  research  it  was  not  expected  that  all  of  these  topics 
could  receive  detailed  attention.  Nevertheless  it  was  hoped  that  a  solid 
base  for  further  work  and  in  particular  cooperation  between  the  Marseille 
and  Arizona  groups  could  be  firmly  established.  This  report  described  the 
significant  progress  that  has  been  made  in  both  of  these  respects.  Of  the 
scientific  objectives,  (a)  and  (b)  were  studied  in  some  detail  but  only  a 
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beginning  could  be  made  in  respect  of  (c),  (d)  and  (e). 


OPTICAL  CONSTANT  MEASUREMENT  AFTER  DEPOSITION 


The  calculation  method  is  described  by  Borgogno  et  al1.  Briefly  it 
consists  of  the  measurement  over  a  wide  wavelength  range  (350-800nm 
usually)  of  the  parameters  T  and  R,  the  transmittance  and  reflectance, 
respectively,  of  the  thin  film  on  its  substrate.  These  parameters  are  then 
used  to  calculate  the  refractive  index,  the  extinction  coefficient,  the 
thickness  and  the  degree  of  inhomogeneity  of  the  layer.  The  model  used  in 
this  determination  exhibits  a  linear  variation  of  refractive  index 
throughout  the  layer  while  the  extinction  coefficient  remains  constant. 
Except  for  the  thickness  all  of  the  parameters  show  dispersion  with 
wavelength. 

Work  under  this  heading  has  been  in  progress  for  some  considerable 
time  in  Marseille.  The  present  contract  allowed  the  completion  of  a  wide 
survey  of  optical  materials  and  the  study  in  greater  detail  of  titanium 
oxide,  one  of  the  most  important  refractory  oxides  for  the  visible  and  near 
infrared  regions.  The  results  are  given  in  Table  I  extracted  from 
Borgogno2. 

An  interesting  study  involved  scandium  oxide,  a  material  that  may  be 
useful  in  the  near  ultraviolet  and  visible  regions  but  about  which  little  Is 
known.  Multiple  determinations  of  the  optical  constants  of  scandium  oxide 
thin  films  that  were  deposited  simultaneously  on  identical  substrates  by 
Optical  Coating  Laboratory  Inc  (the  Marseille  substrates  were  slightly 
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MATERIAL 

A 

A' 

B 

B' 

C 

C‘ 

ZINC  SULPHIDE 

2,2524 

1 ,5754X1 04 

5,4224X10° 

0 

0 

0 

SILICA 

1 ,4625 

3,0693X1 03 

-2,0190X10® 

0 

0 

0 

ALUMINA 

1 .6000 

7,7495X1 03 

-3,9484X1 03 

-0,02 

0 

0 

MAGNESIUM 

1,33619 

8.1846X103 

-5,7723X10® 

FLUORIDE 

0 

0 

0 

ZIRCONIUM 

2,0125 

1.4280X104 

- 1 ,8772X1 07 

OXIDE 

-0,029 

-1,2X10-4 

0 

ZIRCONIUM 

2,0569 

-7.1632X103 

1 ,8979X1 07 

OXIDE 

-0,065 

-1,4X1  O'4 

1 ,9X1 0"7 

ZIRCONIUM 

2,0092 

-5,3327X1 03 

2,0786X10° 

OXIDE 

-0,214 

8,8X1 0-4 

-8,4X1 0-7 

TANTALUM 

1 ,2 1 54X 1 04 

2,2166X10° 

OXIDE 

-8.1X101 

4,2X1 O4 

Table  1:  Values  of  A,  B(nm2)  and  C(nm4),  the  Cauchy  coefficients  in  the 
expansion  for  refractive  index: 

n  -  A  ♦  B/A2  ♦  C/A4 

and  of  A',  B'(nm2)  and  C'(nm4),  the  Cauchy  coefficients  in  the  similar 


expansion  for  An/nm#an 


^inner  ^outar^^ 


mean 
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MATERIAL 

A 

B 

C 

A' 

B' 

C' 

TANTALUM 

2,0841 

3,8096X1 03 

2,2336X1 09 

OXIDE 

-0,017 

5.2X101 

-6,7X1 04 

YTTRIUM 

1 ,7493 

1.1815X104 

2,0971X10® 

OXIDE 

0 

0 

0 

YTTRIUM 

1 ,7525 

1 ,7330X1 04 

-2,9163X10® 

OXIDE 

0,028 

-3,4X1 0-5 

0 

YTTRIUM 

MEM 

1 ,2276X1 04 

-7,0554X1 03 

OXIDE 

0 

0 

YTTRIUM 

1.1223X104 

3,0906X1 07 

OXIDE 

-7,7X10-5 

0 

YTTRIUM 

5,2733X1 03 

4,970 1X1 0s 

OXIDE 

0 

0 

HAFNIUM 

1 ,8998 

1 ,7756X1 04 

-5,7327X10® 

OXIDE 

-0,054 

1,9X10-5 

0 

HAFNIUM 

1 ,9165 

2.1981X104 

-3,2759X10® 

OXIDE 

-0,054 

1,8X10-5 

0 

Table  1  continued 
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MATERIAL 

A 

B 

C 

A' 

B' 

C' 

HAFNIUM 

1 ,8864X1 04 

-2,2999X1 08 

OXIDE 

0 

0 

TITANIUM 

2,0044 

8,7927X1 03 

5,3890X1 09 

OXIDE 

-0,040 

0 

0 

TITANIUM 

2,0357 

3,8484X1 04 

3.7103X109 

OXIDE 

-0,070 

0 

0 

TITANIUM 

2,2658 

-3.4617X104 

1 1.1686X1010 

OXIDE 

-0,040 

0 

0 

TITANIUM 

2,3126 

-3,7938X1 04 

1,241  1X1 010 

OXIDE 

-0,055 

0 

0 

TITANIUM 

- 1 ,41 59X1 04 

9.8190X109 

OXIDE 

-2.69X1 0-4 

0 
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wedged  and  of  different  shape  from  the  others  but  of  the  same  material) 
were  arranged  to  be  made  by  seven  different  groups  including  both 
Marseille  and  Arizona,  who  reported  their  results  at  the  Optical  Society  of 
America  Fall  Meeting  in  New  Orleans  at  the  end  of  October  1983.  A  joint 
paper3  has  been  published  on  the  investigation.  The  Marseille  group  used 
their  normal  techniques  for  reduction  of  the  results:  the  Arizona  group 
made  similar  measurements  of  T  and  R  but  used  a  completely  different 
analytical  approach  that  was  developed  from  a  technique4  for  extracting 
the  optical  constants  of  homogeneous  films  from  measurements  of  T  only. 
The  results  of  the  Marseille  and  Arizona  measurements  agreed  remarkably 
well  as  evidenced  by  Table  2.  The  Arizona  method  uses  the  envelopes  of 
the  transmittance  and  reflectance  characteristics  and  has  been  further 
developed  into  an  in  situ  measurement  technique,  to  be  described  later,  for 
inhomogeneous  film  characterization  using  T  measurements  only. 

An  extensive  study  of  titanium  oxide  has  been  made  in  an  attempt  to 
correlate  the  optical  constants  with  the  conditions  under  which  the  layers 
were  deposited.  This  has  proved  to  be  an  exceedingly  complicated  topic 
and  no  simple  theory  has  so  far  emerged.  The  results  are  extensive  and  are 
still  not  completely  analyzed  and  they  are  specific  to  one  particular 
evaporator.  A  summary  of  the  results  has  been  published  by  Borgogno  et 
al5  and  is  shown  in  figure  1.  The  text  of  this  paper  is  given  as  appendix  1. 

It  is  worth  emphasising  that  optical  inhomogeneity  is  found  in  almost 
all  thin  films  and  this  has  important  implications  for  the  production  of 
high-performance  optical  coatings  especially  since  the  model  used  in  the 
design  and  calculation  of  such  multilayers  is  the  common  parallel-sided 
homogeneous  slab. 
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MEAN  REFRACTIVE  INDEX 


Figure  !.  The  refractive  index  of  titanium  dioxide  layers  depends 
strongly  on  the  deposition  conditions  referred  to  as  a  number  near  each 
curve.  The  values  of  n  as  a  function  of  wavelength  are  deduced  from 
reflectance  and  transmittance  measured  in  air. 
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WAVELENGTH 

OUTER  INDEX 

INNER  INDEX 

nm 

ENVELOPE 

SPECTROMETRIC 

ENVELOPE 

SPECTROMETRIC 

THIN 

333 

1,830 

1,807 

2,038 

2,055 

FILM 

415 

1,802 

1,814 

1,980 

1,960 

551 

1,775 

1,829 

1,914 

1,908 

THICK  315 

1,810 

1,81  1 

2,021 

2,017 

FILM 

344 

1,804 

1,804 

1,999 

2,014 

377 

1,796 

1,790 

2,000 

1,995 

421 

1,787 

1,787 

1,993 

1,975 

478 

1,802 

1,790 

1,957 

1,941 

555 

1,810 

1,797 

1,895 

1,909 

668 

1,819 

1,807 

1,847 

1,875 

Table  2:  Comparison  of  measurements  of  the  refractive  indices  of 
scandium  oxide  made  on  similar  samples  after  deposition  at  Optical 

Sciences  (envelope)  and  Marseille  (spectrometric). 
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In  situ  measurements  are  important  for  several  reasons.  Accurate 
control  of  the  process  demands  an  accurate  knowledge  of  the  optical 
constants  of  the  material  as  it  is  actually  deposited.  These  are  seldom 
exactly  the  same  as  post-deposition  values  which  are  frequently  perturbed 
by  moisture  adsorption.  As  we  have  already  commented,  inhomogeneity  of 
optical  constants  is  the  norm  rather  than  exception  and  there  are  great 
difficulties  in  the  estimation  of  the  actual  profile  of  refractive  index 
through  the  layer  from  post-deposition  results  only.  Any  serious  attempt 
to  relate  microstructure  to  optical  constants  demands  a  knowledge  of  the 
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index  profile  and  hence  more  information  than  can  be  obtained  from 
post-deposition  measurements. 

In  situ  measurements  have  been  an  important  objective  of  the 
Marseille  group  for  at  least  ten  years  and  considerable  progress  has  been 
made  in  the  development  of  equipment  and  especially  in  the  signal 
processing  aspects  of  the  problem.  Measurements  of  film  transmittance 
during  deposition  over  a  wide  spectral  range  are  recorded  using  a  rapid 
scanning  spectrometer.  The  latest  version  of  the  apparatus  uses  concave 
holographic  diffraction  gratings  with  photodiode  arrays  of  the  Reticon 
type  as  detectors  and  the  signals  are  processed  by  PDP1 106  computers. 
Construction  of  an  in  situ  monitoring  system  began  in  Arizona  some  three 
years  ago  under  a  DARPA  contract  through  NWC.  This  system  differs  from 
that  developed  in  Marseille  in  that  a  CCD  detector  array  is  used  rather  than 
the  Reticon  and  the  computer  is  an  IBM  PC6.  An  original  intention  of  the 
present  contract  was  the  provision  of  the  salary  of  a  post-doctoral 
assistant  in  the  Arizona  group  who  would  be  chosen  from  Marseille  so  that 
rapid  profit  could  be  derived  from  the  Marseille  experience.  In  the  event  it 
proved  difficult  to  arrange  for  this  activity  to  be  funded  through  Marseille 
and  a  separate  AFOSR  contract  was  placed  in  the  USA  to  cover  it.  A 
separate  report  covering  only  this  aspect  has  been  submitted  to  AFOSR. 
Nevertheless,  information  is  included  here  for  completeness  and  appendix 
2  summarises  the  design. 

In  situ  measurements  have  therefore  been  made  both  in  Marseille  and 
in  Arizona.  While  the  measurements  were  essentially  the  same,  slightly 
different  techniques  for  the  reduction  of  the  measurements  to  actual  index 
profiles  have  been  used. 

1  1 
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The  method  in  use  in  Marseille  relies  on  the  fact  that  the  admittance 
locus  for  a  single  homogeneous  layer  intersects  the  real  axis  at  the 
turning  values  of  reflectance.  These  points  of  intersection  are  spaced 
apart  by  quarterwave  thicknesses.  Layers  that  exhibit  only  slight 
inhomogeneity  depart  to  a  negligible  extent  from  this  condition.  The 
values  of  the  intercepts  with  the  real  admittance  axis  can  be  readily 
calculated  from  the  value  of  reflectance  at  the  turning  values.  Once  these 
values  have  been  established  it  is  a  straightforward  matter  to  calculate 
the  equivalent  homogeneous  dielectric  layer  between  each  and  this  in  turn 
leads  to  a  curve  of  index  against  thickness  of  the  dielectric  layer  with 
points  every  quarterwave.  The  method  is  described  in  greater  detail  in 
two  publications7,8.  One  of  the  materials  studied  was  Ti02  and  some 

results  are  reproduced  in  figure  2.  Further  measurements  were  made  and 
reported  by  Borgogno  et  al5,  given  as  appendix  I.  These  comprised  two 
sets  of  measurements  and  calculations.  First  the  inhomogeneity  was 
measured  in  the  normal  way  during  deposition.  These  results  are  shown  in 
figure  3  (figure  6  of  appendix  1 ).  Next,  after  deposition,  the  same  layers 
were  removed  from  the  coating  chamber,  R  and  T  measurements  in  air 
were  made,  and  their  optical  constants  calculated,  figure  4  (appendix  1 , 
figure  4).  The  shift  of  index  between  both  of  these  values  is  shown  in 
figure  5  (appendix  1,  figure  7).  The  differences  between  the  two  values  of 
index  can  largely  be  attributed  to  the  entry  of  moisture  in  the  layers  but 
the  details  of  the  difference  curves  indicate  that  more  work  is  necessary 
to  explain  them  fully. 

The  Arizona  method  is  similar  to  that  at  Marseille,  but  relies  on  a 
slightly  different  concept.  The  method  is  described  in  greater  detail  by 
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Fig.  7.  Differences  between  the  values  of  the  refractive  indices 
measured  in  air  i  Fig.  31  and  those  deduced  from  in  situ  measurements 
1  Fig.  51. 
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Fig  i  Measured  transmittance  of  the  titanium  oiide  layer  'Ref. 
h  of  Fig  ii:  1.  measured  in  situ  at  the  end  of  the  deposition:  2.  mea¬ 
sured  m  air  three  days  after  deposition:  3.  calculated  after  Harris  et 

ai.6 

Here,  the  curve  gives  the  result  recerded  two  days  after 
air  admission,  the  spectral  profile  being  well  stabi¬ 
lized. 

Figure  9  gives  the  same  measurements  made  on  a 
sample  'Ret.  2D.  Fig.  3)  the  mean  index  of  which  is  much 
ower.  The  effect  of  air  admission  is  considerable. 

A  first  interpretation  of  these  differences  is  made 
ising  the  calculation  method  developed  by  Harris  et  al  * 
starting  from  the  values  n\\)  determined  in  vacuo,  we 
alculate  the  values  of  the  packing  coefficient  and  of  the 
efractive  index  of  the  bulk  material  that  explain  the 
troperties  measured  in  air.  In  these  calculations,  the 
nterstitial  medium  is  water,  the  index  of  which  nu  \  A) 
5  well  known.  This  leads  to  p  =  0  92  for  the  layer  of 
tef.  6:  the  corresponding  calculated  profiles  are  plotted 
a  Figs.  8  and  9.  Agreement  with  experimental  values 
>  far  from  good. 

Knowing  that  our  technique*  of  index  determination 
tat  take  account  of  mhomogeneity  are  proved  to  be 

70  AOPL.ED  OP'-C 3  cci  23.  Nc  20  IS  Oclouer  1384 


t  r«Q«m:  t  ta  qc« 


Fig.  9.  Measured  transmittance  of  the  titanium  oxide  layer  'Ref. 
20  of  Fig.  3):  1.  measured  in  situ  at  the  end  of  the  deposition:  2. 
measured  in  air.  three  days  after  deposition:  3.  calculated  alter  Harris 

et  aiT 

correct  for  a  quasi-perfect  explanation  of  measured 
optical  properties,  we  hope  to  be  able  to  introduce  this 
additional  parameter  and  so  complete  the  calculation 
methods  of  Macleod  for  a  better  interpretation  of  water 
adsorption  in  layers. 

The  authors  wish  to  thank  the  DRET  which  has 
supported  this  study.  Cooperation  between  Marseille 
and  Tucson  was  especially  enhanced  by  an  Air  Force 
Office  of  Scientific  Research  lAFSOR  84-00241  grant. 

This  paper  was  presented  at  the  third  topical  meeting 
on  Optical  Interference  Coatings.  Monterev,  Calif. 
17-19  Apr.  1984. 
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studied  .ire  —  3 ' 1  nm  thick,  their  extir.cion  c<  efficient 
ij  is  always  smaller  than  2  X  L0~'.  The  evaporated 
material  is  the  oxide  TiO.  We  cannot  give  till  the  details 
or  the  different  evaporation  conditions  indicated  by  a 
number  near  each  curv  e.  Generally  speaking  we  repeat 
Pulker's  results;  the  index  increases  with  substrate 
temperature  which  varies  from  150°C  lexp.  20)  to300°C 
i exp.  no  Many  other  parameters,  such  as  deposition 
rate,  partial  pressure  of  ox;,  gen.  and  state  oft  he  TiO 
material  in  the  crucible,  can  modify  the  values  of  n. 

Concerning  homogeneity  defects,  the  problem  is  even 
more  complicated.  Measurements  lead  to  the  value 
An  n.  We  can  have  either  a  marked  ;r.  homogeneity 
An  n  s  -l  AT  'exp.  20)  or  practically  homogeneous 
layers  with  An  n  s  G  <exp.  3).  The  minus  sign  indicates 
a  decrease  of  refractive  index  from  the  substrate  i Fig. 
2)  m  -  n  -  din  is  systematically  negative  or 

null.  A  general  idea  ut  results  on  the  different  samples 
is  given  .n  Fig.  4.  Around  the  values  of  n.  vertical  lines 
show  the  amplitude  of  n,  -  n  It  seems  that  oxygen 
partial  pre-*ure  and  preparation  >e  g.  cleaning)  con¬ 
dition.-.  of  the  substrate  'see  B-^  and  ~tiica>  have  a 
greater  influence  on  the  inhoniugeneity  ratio  than  does 
temperature. 

8.  In  Situ  Measurements 

The  results  of  in  <itu  measurements  must  be  taken 
with  '(.me  reservations  because  we  have  only  Tar.s- 
mittance  measurements  which  are  less  precise  than 
those  made  in  air.  But,  mainly  tn  iac::u.  during  the 
growth  of  the  layer,  absorption  can  he  disturbing. 
.Sample  1  4  that  gave  the  higher  index  in  air  <see  Fig.  3) 
is  not  correctly  determined,  coefficient  k  being  much 
greater  than  A  X  10“ j.  For  other  samples  it  does  not 
seem  that  absorption  falsities  our  mea-ur..-ments. 

Tlie  me  m  value  of  me  is  tired  in  :  :•  ■■>.  is  given  in 
Fig.  5.  We  also  have  the  law  of  variate  a  of  index  vs 
iaver  thwkn.ess  f> »r  some  thirty  wavelerg’-hs.  For  in 
ev.i;  ■  r.o'mn  with  a  constant  deposition  rare,  the  varia- 
t.  :i  : -v  of  -i i ;  i  is  imear.  As  for  Fig.  -3  representing 
me  i-ur-men's  m  air.  we  represent  the  inhoinoger.viiy 


Fig.  3.  Retractive  index  n  of  titanium  dioxide  layers  of  Fig.  3  cal¬ 
culated  from  in  n:u  measurements  during  deposition. 
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Fig.  6.  Inhomogeneity  An  of  titanium  dioxide  layers  of  Fig.  3  cal¬ 
culated  from  .n  situ  measurements  during  deposition;  the  index 
systematically  decreases  from  the  substrate. 


detected  in  vacuo  by  plotting  the  extreme  values 
<  n.  —  a., )  of  the  samples  i  Fig.  6). 

In  a  general  way.  inhomogeneity  in  vacuo  is  greater 
than  that  observed  in  air.  We  shall  see  that  the  varia¬ 
tions  in  the  mean  index  are  also  significant. 

IV.  Effect  of  Air  Admission  After  Deposition 

First,  let  us  compare  Figs.  3  and  5  to  calculate  the 
shift  observed  in  the  mean  value  of  the  index.  Figure 
7  gives  the  differences  between  indices  vs  wavelength. 
Thi-  iifferer.ee  varies  with  sample  from  0.05  to  0.15. 
We  can  try  to  separate  effects  due  to  the  temperature 
change  'between  the  two  measurements)  and  those  due 
to  the  phenomenon  of  water  adsorption. 

Figure  3  gives  transmittance  vs  wavelength  for  the 
-ample  of  Ref.  6,  in  the  400-1000-nm  range.  This 
measurement  was  made  immediately  following  dspo- 
sition  with  a  substrate  temperature  of  —30t)3C.  The 
return  to  room  temperature  gives  a  slight  modification 
of  ‘he  -pec'ral  profile  that  can  just  he  detected  with  this 
spec To/r.eter.  But  as  -oon  as  air  is  admitted  into  the 
ch  ur.hr:,  the  spectral  profile  changes  spectacularly. 
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'ig.  I.  Model  chosen  to  represent  the  inhomogeneous  layer  between 
ledia  of  indices  n ,  (substrate)  and  n„  lair).  The  index  of  the  layer 
aries  from  n,  at  the  inner  surface  to  n0  at  the  outer.  Inhomogeneity 
i  given  by  An/n  which  is  negative  for  the  index  decreasing  linearly 
from  substrate  to  outer  surface. 

I.  In  Situ  Measurements  During  Deposition  of  the 
aver 

In  spite  of  experimental  difficulties  which  appreciably 
imit  the  precision  of  index  determination,  knowledge 
f  optical  constants  of  layers  in  their  production  con- 
iitions  is  particularly  useful,4  especially  because  these 
alues  are  sometimes  different  from  those  measured  at 
lormai  temperature  and  pressure. 

The  principle  of  the  measurement  is  as  follows:  The 
vaporating  chamber  is  equipped  with  a  rapid  scanning 
pectropnotometer  which  allows  us  to  record  the  evo- 
ltion  of  the  transmittance.  Figure  2  gives  an  illus- 
•ation  of  the  observed  phenomenon.  Calculations  are 
lade  systematically  for  thirty-two  wavelengths  regu- 
rly  distributed  in  the  0.4-1. 1-urn  range.  In  these 
ilculations,  the  extinction  coefficient  is  assumed  to  be 
•gligible.  It  is  then  sufficient  to  extract,  for  each 
avelength.  the  extreme  values  of  the  transmission 
•tected  during  the  deposition. 

Starting  from  two  consecutive  values  T ,  and  TI  +  1, 
at  is,  a  maximum  and  a  minimum,  It  is  possible  to 
construct  the  circle  of  admittance5  and  to  determine 
e  mean  value  of  the  refractive  index  of  the  thickness 
ded  when  the  layer  thickness  changes  from  i A*/4  to 
4  1  )  A*/4.  So,  for  layers  with  large  optical  thicknesses 
veral  A/4),  we  can  determine  the  inhomogeneities, 
ce  we  have  one  value  for  n  for  each  value  of  the  op¬ 
al  thickness  z  equal  tor  *  <i  +  1/2) /4. 

>Ve  shall  examine  some  experimental  results  derived 
m  these  two  measurement  methods  to  emphasize  the 
tematically  occurring  differences  and  explain  their 
?in. 

Inhomogeneity  and  Index:  Dependence  on  the 
mation  Condition 

V'e  have  systematically  studied  the  optical  constants 
itanium  oxide  layers  to  show  how  they  depend  on  the 
osition  conditions.  Many  studies  of  this  topic. 

i 


*00.  KXM.  ^'ihlcknfcssinm) 


Fig.  Z.  Three-dimensional  plot  representing  expected  development 
of  the  transmittance  over  the  400-1000-nm  spectral  region  during  the 
deposition  of  a  layer  of  titanium  oxide  6Aq/4  thick  with  Xo  3  350 

nm. 


Fig.  3.  Refractive  index  of  titanium  oxide  layers  depends  strongly 
on  the  deposition  conditions,  see  the  number  near  each  curve.  The 
values  of  ff  as  a  function  of  wavelength  are  deduced  from  reflectance 
and  transmittance  measured  in  air. 


especially  by  Pulker  and  Ritter,1  have  been  published 
and  they  give  an  idea  of  the  incredible  complexity  of  the 
behavior  of  this  material  evaporated  in  a  reactive  at¬ 
mosphere.  The  results  given  here  concern  refractive 
indices  measured  in  vacuo  and  these  same  indices 
measured  after  cooling  and  exposure  to  atmospheric 
humidity.  They  demonstrate  the  complexity  of  the 
problem  rather  than  leading  to  a  better  understanding 
of  it. 

Nevertheless,  we  note  that  knowledge  of  these  values, 
along  with  the  control  of  the  layer  deposition  conditions, 
is  absolutely  necessary  to  ensure  the  accuraie  manu¬ 
facture  of  multilayer  stacks.  We  can  give  only  a  general 
idea  of  our  experimental  results. 

A.  Measurements  in  Air 

Figure  3  shows  the  dispersion  law  vs  wavelength  of 
the  mean  value  of  index  n,  the  layers  being  measured 
in  air. 

We  have  chosen  a  sample  group  of  layers  to  show  the 
disparity  of  the  values  normally  used.  All  the  layers 
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Tae  tact  that  the  uptual  vharacteristici  •  r  'h.p.-ftim  .T.aterials  are  generally  different  from  those  of  the  same 
materials  in  hulk  turm  .?  veil  known.  The  -inferences  depend  very  much  on  the  conditions  in  which  the  de¬ 
position  has  been  carried  .'Ut.  A  gad  .."irrscanding  of  these  differences,  their  causes,  and  the  influence 
ot  deposition  parameters  :.s  vital  it  we  are  to  be  able  to  improve  coating  quality.  We  have  developed  two 
complementary  men's  d»  with  the  objective  of  deriving  information  on  the  index  of  refraction  and  its  varia¬ 
tion  throughout  the  tnii.kn.ess  ot  the  film.  Perceptible  optical  inhomogeneity  is  normally  present  and  appre¬ 
ciable  inhomogeneity  is  frequently  present  in  thin  films.  Such  inhomogeneity  is  usually  associated  with 
layer  micro-structure.  The  tirst  is  a  p'.scdeposrion  technique  that  makes  use  of  measurements  in  air  of  the 
transmittance  and  reflectance  of  the  layer  under  study  over  a  '.vide  wavelength  region.  The  second,  in  con¬ 
trast,  makes  use  of  ;-■»  >;.‘u  measurements,  that  is  measurements  made  under  vacuum  and  during  the  actual 
deposition  ot  the  layer.  We  shall  show  with  the  help  of  several  examples  that  the  two  methods  lead  to  results 
that  are  consistent  and  demonstrate  the  existence  in  deposited  materials  of  an  inherent  variation  of  the 
index  of  retraction  normal  to  the  surface.  The  thermal  sensitivity  of  the  layer  properties  and  their  tendency 


to  adsorb  atmospheric  moisture  must 
techniques  can  be  explained. 


I.  Introduction 

In  the  literature1  are  many  papers  on  the  influence 
of  evaporation  conditions  on  the  stoichiometry  and 
microstrucrure  of  evaporated  iayers.  For  example,  for 
numerous  materials  the  complex  index  of  refraction 
obviously  depends  on  these  conditions  and  it  is  impor¬ 
tant  to  take  account  of  this  fact  in  solving  the  problems 
of  manufacturing  optical  muitidielectric  niters..  In  fact, 
in  most  cases,  the  values  of  refractive  indices  are  de¬ 
duced  from  accurate  spectrophotometric  measurements 
made  on  the  layers  after  they  are  taken  out  of  the 
evaporating  chamber. : 

An  accurate  control  of  the  growth  of  the  layers 
involves  in  ,ifa  knowledge  of  indices,  that  is.  as  they  are 
actually  growing.  These  values  can  be  different  from 
those  observed  in  air,  at  normal  pressure  and  temper¬ 
ature. 

We  shall  briefly  describe  the  techniques  used  for 
determining  the  indices  from  the  different  measure- 
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be  taken  into  account  before  the  residual  differences  between  the  two 


ments  and,  using  some  examples,  we  shall  show  how 
complicated  are  the  experimental  results  given  by  this 
study. 

II.  Techniques  for  Determining  Optical  Constants 

We  can  distinguish  two  types  of  method:  measure¬ 
ments  made  in  air  and  those  made  in  vacuo  during  the 
deposition  of  the  layers. 

A.  Measurements  in  Air3 

The  principle  of  the  method  of  determination  of  op¬ 
tical  constants  is  classical:  reflection  and  transmission 
factors  of  a  layer,  deposited  on  a  prismatic  substrate  of 
well-known  index,  are  measured.  Our  measurements 
are  made  over  a  wide  spectral  range:  from  0.4  and  1.1 
urn,  In  most  cases  the  hypothesis  of  a  homogeneous 
layer  is  not  sufficient  to  correctly  explain  the  measured 
optical  properties.  For  the  calculation,  the  layer  is 
therefore  assumed  to  have  a  constant  index  gradient  in 
its  thickness  and  the  three  following  parameters  are 
determined  conjointly: 

the  mean  index  of  complex  refraction  n‘  =  n  —  ik, 
which  is  a  function  of  the  wavelength 

the  relative  variation  of  index  in  the  thickness  of  the 
laver  An,  ft  with  the  sign  convention  given  in  Fig.  1; 

the  geometrical  thickness  d. 

Measurements  and  calculations  are  made  automatically, 
and  we  have  done  a  complete  and  systematic  study  in 
our  laboratory  to  establish  the  dependence  of  n.  k,  and 
An  n  on  the  evaporating  conditions  of  materials  gen¬ 
erally  u.zed  in  the  visible  spectral  range. 
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FIGURE  CAPTIONS 

1 .  The  refractive  index  of  titanium  dioxide  layers  depends  strongly  on 
the  deposition  conditions  referred  to  as  a  number  near  each  curve. 

The  values  of  n  as  a  function  of  wavelength  are  deduced  from 
reflectance  and  transmittance  measured  in  air. 

2.  Titanium  oxide.  Values  of  refractive  index  as  a  function  of  thickness. 
The  results  are  derived  from  the  measurements  recorded  during  the 
deposition  for  the  wavelengths: 

479.0nm  ▲  555.8nm  □  594.2nm  ■  632.6nm  o  728.6nm 
t  824.6nm 

The  optical  thicknesses  4  \/4  and  8  ^*/4  are  located  on  the  thickness 
axis.  -  5  So  nm 

3.  The  inhomogeneity  An  of  titanium  dioxide  layers  of  fig.  I  calculated 
from  in  situ  measurements  during  deposition;  the  index 
systematically  decreases  from  the  substrate. 

4  Homogeneity  defects  of  the  layers  referred  to  in  figure  1 . 

An”nmsi<j«  ‘  noutside  <s  Plotted  as  3  function  of  wavelength;  the  index 
systematically  decreases  from  the  substrate. 

5  Difference  between  the  mean  values  of  the  refractive  indices 
measured  in  air  and  in  situ. 
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Monitoring  of  optical  thin  films  during  deposition 
Proc  SPIE,  Vol  40 1 ,  pp  74  -  82,  1 983 

9.  F  Flory,  B  Schmitt,  E  Pelletier  and  H  A  Macleod 

Interpretation  of  wide  band  scans  of  growing  optical  thin  films  in 

terms  of  layer  microstructure 

Proc  SPIE,  Vol  401,  pp  109  -  116,  1983 
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CONCLUSIONS 

Techniques  for  the  measurement  of  thin-film  optical  constants  both 
during  and  after  deposition  have  been  developed  at  the  Ecole  Nationale 
Superieure  de  Physique,  Marseille,  and  the  Optical  Sciences  Center, 
Arizona.  These  techniques  have  been  used  to  study  the  inhomogeneity  of 
films,  especially  titanium  dioxide,  and  the  way  in  which  the  optical 
constants  change  with  moisture  adsorption  on  exposure  to  the  atmosphere. 
These  aspects  of  film  behaviour  are  directly  related  to  film 
microstructure.  The  collaboration  between  the  two  groups  has  been 
fruitful,  particularly  in  respect  of  the  contribution  that  was  made  to  the 
work  at  Optical  Sciences  by  B  Bovard  from  Marseille. 
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Bovard  et  al9  which  is  given  as  appendix  3.  The  layer  is  assumed  to  have 
slight  inhomogeneity  that  varies  smoothly  throughout  it.  In  the  case  of  a 
lossless  dielectric  film  this  leads  to  a  straightforward  and  slight 
modification  of  the  normal  matrix  expression  for  the  film.  Small  losses 
can  be  included  by  altering  the  phase  thickness  to  include  a  complex 
refractive  index.  An  envelope  method4  similar  to  that  employed  in  the 
post-deposition  determination  of  film  indices3  is  then  used  to  extract  the 
inner  and  outer  values  of  refractive  index  at  each  stage  of  film  growth. 
The  principal  problem  is  in  drawing  the  envelopes.  Numerical  techniques 
for  smoothing  the  curve  of  transmittance  against  thickness  are  first 
employed  and  then  parabolic  expressions  are  used  to  generate  the  envelope 
curves.  Measurements  on  a  typical  titanium  dioxide  layer  prepared  under 
good  conditions  give  for  a  total  geometrical  thickness  of  670nm,  an  inner 
refractive  index,  that  is  the  value  next  to  the  substrate,  of  2. 135  and  an 
outer  refractive  index  of  1.794.  The  variation  of  both  these  parameters 
with  wavelength  is  given  in  appendix  3,  figure  4. 

The  method  implies  the  assumption  that  the  layer  is  sufficiently 
stable  during  the  measurement  that  the  values  of  the  part  already 
deposited  do  not  change.  This  can  in  fact  be  used  as  a  test  for  stability. 
Other  measurements  of  the  index  profile  of  a  layer  prepared  under  poor 
oxidising  conditions,  appendix  3,  figure  5,  show  considerable 
discontinuities  in  extinction  coefficient  that  even  becomes  negative.  This 
indicates  a  lack  of  stability  of  such  a  layer. 

This  work  represents  simply  a  beginning.  We  must  carry  out  further 
comparisons  of  the  two  methods.  This  can  even  be  done  on  the  same  data. 


Wave 1 ength (nm) 

ITTTANIIIM  OXIDE  (  measurements  in  air  )  | 


e  4.  Homogeneity  defects  of  the  layers  referred  to  in  figure  1 . 
An.n  -  n  „  is  plotted  as  a  function  of  wavelength;  the  index 

An  outstd* 

smatically  decreases  from  the  substrate. 
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Figure  3.  The  inhomogeneity  An  of  titanium  dioxide  layers  of  fig.  1 
calculated  from  in  situ  measurements  during  deposition;  the  index 
systematically  decreases  from  the  substrate. 
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Figure  2.  Titanium  oxide.  Values  of  refractive  index  as  a  function  of 
thickness.  The  results  are  derived  from  the  measurements  recorded  during 

the  deposition  for  the  wavelengths: 

?  479. Onm  ▲  555.8nm  □  594.2nm  ■  632.6nm  o  728.6nm 

•  824.6nm 

The  optical  thicknesses  4  A./ 4  and  8  X/4  are  located  on  the  thickness 
axis.  (  -  SSc  f'rT>>) 
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ABSTRACT 

A  scanning  monchromator  system  -for  the  monitoring  of  thin  film 
deposition  in  a  box  coater  is  described.  The  system  employs  data 
from  both  a  quartz  crystal  oscillator  and  a  wide  band  transmission 
spectrometer.  The  spectrometer  uses  a  holographic  grating  as  its 
dispersive  element  and  a  CCD  array  to  collect  the  data.  All  data  is 
sent  to  a  microcomputer  where  the  information  is  displayed,  stored, 
and  analyzed.  Several  applications,  including  measurement  of 
optical  constants  of  inhomogeneous  films  and  characterization  of 
moisture  adsorption,  are  discussed. 
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1.  INTRODUCTION 

After  an  optical  filter  has  been  satisfactorily  designed,  it 
must  be  implemented  in  a  production  facility.  For  anything  but  the 
simplest  of  multilayer  stacks,  this  involves  the  selection  of  the 
proper  process  parameters  for  a  specific  set  of  materials  in  a 
particular  vacuum  chamber.  The  length  of  time  required  to  control 
this  process  fully  depends  upon  the  efficiency  and  fidelity  of  the 
monitoring  techniques  available  to  the  operator.  In  general,  there 
are  two  means  of  monitoring  film  deposition: 

1)  Optical  monitoring:  This  technique  is  more  appropriate  in 
realizing  coatings  consisting  of  quarterwave  layers  by  detecting  the 
extrema  of  transmission  or  reflectance  at  a  particular  wavelength. 
This  method  is  extremely  stable  for  the  control  wavelength  and  only 
slightly  less  stable  around  it.  1  Therefore,  it  is  used  for  coat¬ 
ings  designed  for  performance  over  a  narrow  wavelength  region. 

2)  Physical  Mass  Monitoring:  By  observing  the  natural  reso¬ 
nance  frequency  of  a  quartz  crystal,  it  is  possible  to  determine  the 
mass  deposited  on  its  surface;  assuming  densities  for  the  film 
materials,  one  can  compute  the  thickness  of  the  layers.  With  this 
technique  one  can  monitor  layers  of  any  thickness  with  high  sen¬ 
sitivity.  Unfortunately,  since  it  does  not  include  any  measurements 
of  the  optical  performance,  it  does  not  provide  the  stability  of 
optical  monitoring,  and  demands  accurate  calibration  and  reliable 


density  data. 
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In  this  paper,  we  will  describe  a  scanning  monochromator 
system  which  employs  both  of  these  techniques  in  parallel  with  the 
added  advantage  of  measuring  the  optical  transmission  of  the  sample 
over  a  wide  wavelength  range.  We  continue  by  discussing  relevant 
•features  of  the  system  and  considering  some  of  its  applications. 

2.  IMPLEMENTATION 

It  is  appropriate  to  begin  by  mentioning  that  our  system  was 
inspired  by  one  built  by  the  group  led  by  Pelletier  in  Marseilles, 

2  q 

France.  *  Our  scanning  monochromator  system  was  intended  to 
augment  the  capabilities  of  a  Balzers  760  box  coater,  which  was 
delivered  with  an  automated  process  controller  (Balzers  Model  KB 
101)  based  on  a  quartz  crystal  monitor  and  a  second,  unautomated, 
single  wavelength  optical  monitor  (Balzers  Model  GSM  210).  The 
front  end  of  the  first  subsystem  was  left  intact;  the  second  sub¬ 
system  was  replaced  with  our  wideband  optical  monitoring  portion  of 
the  scanning  monochromator  system. 

Figure  1  is  a  descriptive  flow  diagram  of  the  scanning  mono¬ 
chromator  system.  To  ensure  that  the  signal  reaching  the  CCD  array 
is  adequate,  the  original  Balzers  light  source  was  replaced  with  a 
much  brighter,  590  W  tungsten-halogen  lamp.  This  source  met  the 
required  criteria:  enough  light  to  saturate  the  detectors  and  suf¬ 
ficient  stability  to  rely  on  a  single  109X  reference  line  for  an 
entire  coating  run.  Unfortunately,  its  spectral  profile,  shown  in 
Figure  2A,  varies  excessively  and  does  not  extend  into  the  lmportan* 
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ultraviolet  region.  To  solve  this  problem,  we  are  currently  repla¬ 
cing  the  tungsten-halogen  with  a  xenon  arc  lamp  which  will  have  both 
a  flatter  spectrum  and  output  in  the  ultraviolet,  shown  in  Figure 
2B. 

Immediately  adjacent  to  the  source,  the  beam  is  modulated  bv  a 
four  sector  chopper,  which  also  provides  a  reference  signal  to  the 
lock-in  amplifier,  part  of  the  detector  electronics  beyond  the  de¬ 
tector.  After  the  chopper,  the  light  is  projected  by  a  lens  into 
the  chamber  via  a  port  in  the  baseplate,  through  the  witness,  or 
reference  sample,  and  back  out  through  a  port  in  the  chamber  roof. 
For  some  experiments,  a  rotating  fixture  moves  more  than  one  sample 
through  the  monitor  beam,  permitting  in  situ  coating  comparisons. 
Figure  3  depicts  the  overall  arrangement  of  the  scanning  monochroma¬ 
tor  with  respect  to  the  original  box  coater. 

After  exiting  the  chamber,  the  beam  is  turned  by  a  flat  mirror 
into  the  mam  scanning  monochromator  optics.  First,  it  is  refocused 
bv  a  lens  onto  a  slit,  which  is  1  cm  high  and  approximately  30  Mm 
wide.  After  the  slit,  the  beam  encounters  the  dispersive  element  of 
the  scanning  monochromator,  a  Jobin-Yvon  holographic,  reflective 
grating,  ruled  at  380  1 ines/mm  and  designed  to  disperse  light  into  a 
spectrum  ranging  from  400  to  308  nm,  of  which  we  use  440  to  300  nm 
due  to  our  light  source.  The  grating  has  sufficient  optical  power 
to  image  the  slit  onto  the  Fairchild  122  CCD  ar^ay,  which  is  reached 
after  a  final  bounce  from  a  flat  folding  mirror.  The  CCD  intercepts 


the  beam  at  an  image  of  the  grating  s  negative  first  order;  tne  one- 
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inch  long  CCD  array  matches  the  one-inch  flat  -field  o-f  the  CCD  at 
that  point.  4  a  view  o-f  this  part  of  the  monochromator  appears  in 
Figure  4. 

The  CCD  arrav  consists  ot  1728  elements;  our  signal  processing 
electronics  averages  sets  ot  ten  adjacent  elements,  providing  us 
with  173  data  points.  These  data  levels  are  sent  on  to  a  dedicated 
IBM-PC,  which  records  the  data  on  tive  inch  floppy  disks  and  dis¬ 
plays  them  on  an  Amdek  video  monitor  for  real  time  feedback  to  the 
plant  operator.  At  the  same  time,  information  from  the  original 
process  controller,  based  on  the  quartz  crystal  monitor,  is  sent 
through  an  interface  module  to  the  I8M-PC. 

A  flowchart  showing  how  the  computer  handies  the  data  appears 
in  Figure  5.  The  IBM-PC  incorporates  a  Tecmar  A/D  board,  which 
accepts  12-bit  data  from  both  the  quartz  crystal  monitor  and  the  CCD 
array.  Although  the  electronics  are  capable  of  running  at  a  rate  of 
four  spectra  per  second,  we  generally  take  one  spectrum  every  three 
seconds.  This  is  quite  adequate,  since  at  our  typical  deposition 
rates  (primarily  of  oxides),  we  deposit  ten  to  twenty  Angstroms  of 
material  in  three  seconds.  The  potential  for  data  rates  at  least  an 
order  of  magniture  higher  would  permit  us  to  monitor  extremely  raDid 
changes,  should  the  need  arise.  Table  1,  below,  summarizes  impor¬ 
tant  characteristics  of  the  scanning  monochromator  monitor. 

Before  each  run,  wavelengths  are  calibrated  bv  inserting  a 
piece  of  didymium  glass  in  the  light  path  and  then  setting  zero 
readings  in  the  data  to  the  didymium' s  Known  absorotion  lines.  For 
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readings  other  than  those  given  by  the  absorption  lines  the  wave¬ 
length  is  determined  by  linear  interpolation  between  known  points. 

We  estimate  an  accuracy  of  2  to  4  nm  over  the  range  we  have  tested 
with  available  spectral  line  sources. 

3.  APPLICATIONS 

Several  examples  of  applications  for  the  scanning  monochromator 
monitor  follow: 

1)  A  sequence  of  transmission  spectra  for  each  run  can  be  stored 
for  later  analysis  of  the  effects  of  various  process  parameters.  A 
360K  double-density  5.25  inch  disk  can  accommodate  about  50  minutes 
of  continuously  monitored  spectral  data. 

2)  The  wideband  transmission  spectra  that  appear  on  the  monitor 
provide  the  plant  operator  with  a  much  broader  view  of  coating  pro¬ 
gress;  should  problems  arise,  the  operators  can  base  their  decisions 
during  deposition  on  a  larger  data  base. 

3)  Monitoring  done  in  situ  permits  testing  and  observation  of 
coatings  without  removing  them  from  the  system.  One  example  of  this 
feature's  utility  appears  in  Figure  6:  It  examines  the  effect  of 
water  adsorption  on  a  filter.  ® 

4)  The  larger  data  base  available  to  the  computer  permits  better 
characterization  of  a  film's  optical  constants.  Figure  7  shows  one 
such  calculation  for  a  film  of  TiO^  as  it  grows.  We  note  that  this 
illustrates  the  dispersion  of  n  and  K  with  thickness  at  a  particular 
wavelengtn.  The  system  can  also  derive  the  dispersion  with  wave- 
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length  at  a  particular  thickness.  ^ 
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TABLE  1:  SYSTEM  PERFORMANCE 


Wavelength  Range 


Wavelength  Resolution 


Resolving  Power 


Etendu  of  Svstem 


449  -  300  nm 


2  nm 


300 

1  x  10  cm2-sr 


Minimum  Transmission 


Z/. 


Signal  Level  o 4  CCD 


7BV.  ot  saturation  at  650  nm 


FIGURE  CAPTIONS 


Scanning  Monochromator  Flow  Diagram 

Spectral  Profile  of  Tungsten  Halogen  Lamp 

Spectral  Profile  of  Xenon  Arc  Lamp 

Appearance  of  Scanning  Monochromator  system 

Top  view  of  Scanning  Monochromator 

Flowchart  of  Computer  Data  Handling  Program 

Example  of  Water  Adsorption  in  a  HO,,  Si02  Fabry-Perot  Filter 

Variation  of  optical  constants  of  Ti02  as  the  film  is  being 

deposited.  Upper  curve  corresponds  to  n,  lower  curve  to  k. 


Scanning  Monochromator  Flow  Diagram. 


IH.  Data  Processing  and  Simulation 


The  envelope  computation  requires  an  accurate  knowledge  of  the  extrema 
together  with  a  method  to  fit  curves  to  these  points.  It  was  thought 
possible  that  the  daca  processing  used  for  the  calculation  of  the 
envelopes  might  introduce  a  bias  in  the  results  given  by  the  program  of 
optical  constant  determination.  However,  we  are  able  to  show  that  it  is 
not  the  case  so  that  our  results  are  free  from  such  effects. 

in.  1.  Smoothing  the  curve  of  transmission  versus  time 
Many  methods  are  available  for  smoothing,  but  our  major  concern  is  to 
3000th  without  distorting  the  curve.  For  this  reason,  we  decided  to  use 
a  filtering  method  based  on  a  finite  impulse  response  filter  calculated 
to  have  a  linear  phase  and  an  extremely  flat  low-pass  band.  We  thus 
avoid  any  distortion  introduced  by  a  nonlinear  phase  and  the  attenuation 
due  to  a  nonflat  passband.  The  filter  does  still  introduce  a  delay 
equal  to  the  derivative  of  its  phase  versus  frequency,  but  this  delay  is 
constant  and  is  corrected  in  the  computation  by  an  entire  shift  of  the 
curve.  Details  of  this  type  of  approach  are  given  in  (9)  which  inspired 
the  present  design  and  so  we  limit  our  description  to  demonstrating  its 
application  to  a  real  signal  in  Figs.  1  and  2. 

III.  2.  envelope  Computation 

The  equations  of  the  envelopes  are  computed  by  segments.  The  first 
group  of  three  points  gives  a  parabola  which  is  used  to  describe  the 
segment  of  the  curve  between  the  first  two  points.  Next  the  first  point 
is  discarded  and  replaced  by  the  fourth  point.  A  new  parabola  is  now 


Provided  we  koov  the  envelopes  of  a  curve  of  transmission,  using  these 
analytical  expressions  we  are  able  to  derive  the  optical  constants.  But 
the  necessity  of  knowing  the  envelopes  as  functions  of  thickness  leads 
us  to  a  very  basic  assumption  concerning  the  stability  of  growth  of  the 
layer.  The  envelopes  calculated  using  the  measured  transmission  versus 
time  are  usable  if  they  do  not  shift  during  the  deposition  process,  which 
means  that  an  earlier  part  of  the  layer  is  not  modified  during  the 
deposition  of  a  later  part. 

It  is  easy  to  foresee  the  importance  of  having  an  accurate  way  of 
determining  extrema  of  the  transmission  curves.  Their  precise 
estimation,  both  value  and  position,  demanded  that  we  consider  some 
numerical  treatments  of  the  raw  data  before  calculating  the  envelopes 
fitting  these  points.  We  present  in  the  next  section  the  data  processing 
techniques  developed  together  with  a  justification  of  their  validity  and 
of  the  expressions  presented  above. 


d 

since  calculation  of  4-  n(z)dz  is  impossible,  we  make  an 

d  Jo 


approximation  and  write  that  °  “  4  n(u)du  where  t  is  the  instant 

io 


when  the  extremum  occurs. 


This  is  equivalent  to  assuming  that  the  rate  of  deposition  (change  in 
thickness  per  unit  time)  is  constant.  The  assumption  is  reasonable 
because  our  deposition  rate  is  automatically  controlled. 


With  regard  to  the  profile  of  extinction  coefficient  we  note  first  that  a 
mean  value  of  the  extinction  coefficient  can  be  derived  each  time  we 
reach  a  new  quarterwave: 


k  -  - 


C  j[  l-(Toax/Tmin)l/*J 


t— j  Log  a(d)  with  a(d)  -  - - 

4»d  C2[  i-KTmax/Tmin)1/* 


*  exp  (-At  ~r~) 


but  since  the  function  a  is  available  at  any  instant  during  the  growth  of 
the  layer,  we  can  also  calculate  its  derivative  versus  time  and  obtain 
an  expression  giving  the  profile  of  extinction  coefficient. 

k(z)  •  ~~  7 — r  ~  where  2  is  the  thickness  of  the  layer. 

At  a(z)  dz 


Because  of  the  derivative,  this  expression  presents  some  sensitivity  to 
errors  but  can  still  be  used  to  obtain  an  indication  of  the  absorption  of 


the  material 


refractive  index,  the  geometrical  thicknesses  and  the  extinction 
coefficient  profile. 


The  expression  of  the  innermost  refractive  index  is  obtained  assuming  n^Q 
-  ani,r  and  d  ■  0  so  that: 


°in  “  where  H- — ~  +  2n0ns 


The  expression  of  the  outermost  refractive  index  is  then: 

„  .  2ninnsao  Tmax-Tmin  +  n  „* . ,  Tmax-Tmin  N, ,„!»/* 

ouC  nijj 2-ng  1  Tmax  Tmin  °|_  10  Tmax  Tmin  ^  ^“in^s 

which  is  calculable  only  if  we  know  the  innermost  index  value. 


Provided  we  assume  the  innermost  index  is  stable  we  can  calculate  the 
profile  of  the  refractive  index  of  an  inhomogeneous  layer.  Note  that  it 
requires  the  knowledge  of  the  substrate  refractive  index,  which  can  be 
measured  independently,  but  it  does  not  require  the  value  of  wavelength. 


We  also  require  the  geometrical  thickness.  Since  we  are  dealing  with 
low-  absorption  materials,  the  extrema  of  transmission  occur  when  the 
optical  thickness  of  the  layer  is  a  multiple  of  a  quarterwave. 

If  m  is  the  order  of  the  extremum  (m*l  indicating  the  first  minimum)  the 
corresponding  geometrical  thickness  is 

rd 

d  *  m  —  where  n  ■  -r  n(z)dz 

.r  *  Jo 


T  -  lfagOp  Qjj  Bout  e25» 

C,2  +  C,2  e46*  +  2ClC,e26*  cos26 
where  C1-(nout+noXng-niin)  C,-(n 
6,*2*nd/ X 


i 


out 


-noXDg-tiiu) 
— 2ikd/i 


In  the  following  we  shall  assume  chat  we  are  dealing  with  a  high  index 
layer  so  that  C,  is  negative..  Then  the  two  expressions  of  the  envelopes 
become: 


16n0nsnl000ute2g« 
max  (Cj  +  C,  e2i»)J 

16no°salnDoute2g» 

010  (Cl  -C,  •2*.)« 

These  two  equations  are  used  to  determine  the  outermost  index  for  each 
instant  during  film  deposition.  They  are  also  used  to  calculate  the 
extinction  coefficient  at  the  quarterwave  points. 

However,  all  these  derivations  are  possible  only  if  we  make  a  very  basic 
assumption  about  Che  stability  of  the  layer  during  deposition.  This 
assumption  is  simply  that  the  innermost  index  does  not  vary  during  the 
growth  of  the  film.  The  profiles  of  the  optical  constants  can  then  be 
derived  and  considered  as  functions  of  thickness  instead  of  time. 

We  shall  not  go  through  the  details  of  the  derivation  and  we  shall  only 
give  the  analytical  expressions  for  the  innermost  and  outermost 


A  non-absorbing  Inhomogeneous  layer*  provided  it  is  reasonably  chick  and 
there  is  appreciable  index  contrast  at  its  boundaries,  can  be  represented 
by  the  characteristic  matrix®1 

"(nin/nout)1/1  cos  5  i  sin  «/(nout  n^)1/*  1 
*<°in  °out)1/*  8in  6  (nout^°in^1/*  C08 

where  nout  is  the  index  at  the  outer  surface  of  the  film 

n^g  is  the  index  at  the  inner  surface  of  the  film,  and 
6  is  the  phase  thickness. 

If  we  assume  that  the  absorption  is  very  small,  the  same  matrix  can  be 
used  with  the  extinction  coefficient  included  in  5  only  so  that 
6  *  2r(n  -  ix)d/l 

r  d 

where  n  is  the  mean  index  of  the  film  (i.e.  n  ■  ~r  n(z)dz), 

11  Jo 
1  f  ^ 

k  is  the  mean  extinction  coefficient  (i.e.  k  ■  -r  \  k(z)dz), 

d  Jo 

X  is  the  wavelength  and 
d  the  thickness  of  the  layer. 

Let  us  assume  that  the  incident  medium  has  a  refractive  index  n0  and 
that  the  substrate  is  nonabsorbant  of  index  ng.  After  some  tedious 
calculation,  we  find  the  expression  for  the  transmission  of  the  coated 


substrate: 


In  Che  case  of  negligible  extinction  coefficients  Che  derlvacion  of  Che 
Tefracdve  index  can  be  achieved  using  an  admittance  circle  method  based 
on  transmission  measurements  during  the  growth  of  the  layer^.  But  since 
our  lncerescs  lie  more  in  slightly  absorbant  materials  such  as  oxides, 
we  have  developed  a  technique  inspired  by  the  envelope  method  described 
by  Manifacier  et  al  (6).  We  have  used  an  inhomogeneous  model  of  a  Chin 
film  to  derive  che  refractive  index  and  Che  extinction  coefficient 
profiles. 

II.  1.  Envelope  method  in  the  Inhomogeneous  case 

Manifacier  et  al^  have  fully  described  the  envelope  method  in  the 
limiting  case  of  an  homogeneous  thin  film  where  only  transmission 
measurements  are  required  for  the  derivation  of  n  and  k.  A 
generalization  to  an  inhomogeneous  model  was  presented  by  Arndt  et  al7 
to  derive  the  optical  constants  from  measurements  of  reflectance  and 
transmittance.  In  these  studies  reflection  and  transmission  are 
considered  as  functions  of  wavelength  and  their  envelopes  are  used  to 
calculate  the  optical  constants.  To  obtain  the  profiles  of  the  optical 
constants,  as  in  this  study,  we  considered  the  envelopes  of  the  curve  of 
transmittance  as  functions  of  thickness  for  a  chosen  wavelength.  We 
were  able  to  obtain  the  dispersion  by  applying  the  same  method  to  a 
range  of  wavelengths. 

We  now  demonstrate  the  principle  of  the  method  and  give  the  analytical 
expressions  for  the  optical  constants. 


L  Introduction 


The  refractive  index  and  the  extinction  coefficient  of  a  thin  film  depend 
upon  the  conditions  of  deposition  and  as  a  consequence  upon  the 
structure  of  the  film  itself.  In  the  case  of  oxides,  inhomogeneities  are 
largely  due  to  the  film's  columnar  structure*  and  to  the  variations  in 
degree  of  oxidation  throughout  the  layer.  The  derivation  of  their 
profiles  as  a  function  of  thickness  is  difficult  once  the  layer  has  been 
deposited.  Furthermore,  most  techniques  developed  to  measure  the 
refractive  index  are  carried  out  under  atmospheric  conditions.  Exposure 
of  a  film  to  the  air  modifies  its  optical  properties:  the  voids  existing 
in  its  structure  tend  to  adsorb  moisture  and  an  oxidation  process  may 
occur  for  a  suboxidized  layer,  changing  the  refractive  index  and  the 
extinction  coefficient. 2. 3  Therefore  a  technique  taking  into  account  the 
evolution  of  the  transmission  of  a  layer  growing  in  vacuo  has  great 
advantages.  To  achieve  such  measurements  we  have  used  a  scanning 
monochromator  system^  which  provides  us  with  the  transmission  over  the 
visible  spectrum  versus  time,  during  deposition.  Using  these  values  we 
have  developed  a  technique  for  deriving  the  profiles  of  the  optical 
constants.  After  a  verification  of  the  optical  constants  determination 
technique  by  computer  simulation,  we  have  applied  the  method  to  various 
layers  of  titanium  dioxide.  This  technique  can  then  be  used  as  a  means 
of  monitoring  the  eifect  on  the  optical  constants  of  a  change  in  any  of 
the  parameters  used  for  the  coating  process. 

11.  Determination  of  the  profile  of  the  optical  constants  of  an 


inhomogeneous  film 
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ABSTRACT 

The  optical  constants  of  a  thin  film  depend  upon  the  structure  of  the 
film  itself.  A  technique,  based  on  transmission  measurements  carried 
out  in  vacuo,  has  been  developed  to  derive  the  profiles  of  the  refractive 
index  and  extinction  coefficient.  The  interpretation  of  the  profiles 
gives  information  on  the  layer  structure  in  vacuo.  The  technique  can  be 
used  as  a  means  of  monitoring  the  variations  of  the  optical  constants 
with  changes  in  the  deposition  parameters.  This  paper  presents  the 
technique,  which  is  based  on  an  envelope  method,  and  gives  some 


experimental  results. 
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Variation  of  optical  Constants  of  Ti02  film  as  the  film  is  being  deposited 
Upper  curve  corresponds  to  n,  lower  curve  to  k. 


SHIFT  IN  WAVELENGTH 
WITH  PUMPDOWN  FOR 
A  NARROW  BAND  FILTER 


Example  of  Water  Adsorption  In  a  TiOj,  SiOj  Fabry -Perot  Filter 


Flowchart  of  Computer  Data 
Handling  Program. 


Top  view  of  Scanning  Monochromator 
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calculated  and  used  to  describe  the  segment  between  the  second  and  third 

points.  This  operation  is  repeated  until  the  entire  curve  is  produced.  ►  -• 
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m.  3.  Verification  and  precision  of  the  technique 

The  validity  of  the  expressions  for  the  refractive  indices  and  the 
extinction  coefficient  has  been  checked  by  some  simulations,  the  results 
of  which  are  presented  here. 

A  program  able  to  handle  variable  optical  constants  has  been  developed 
to  compute  the  evolution  of  the  transmission  of  a  coated  substrate 
during  growth  of  a  thin  film.  The  profiles  of  the  optical  constants  are 
fed  into  the  program  which  divides  the  inhomogeneous  layer  into  a  stack 
of  homogeneous  sublayers  each  of  thickness  Inm.  The  transmittance  is 
then  recorded  on  disk  files  exactly  as  when  the  scanning  monochromator 
is  used.  The  optical  constant  determination  program  is  then  used  to 
compute  the  index  and  estimation  coefficient  profiles.  Their  comparison 
with  the  original  ones  verifies  the  accuracy  of  the  technique.  Our  goal 
was  also  to  show  it  is  possible  to  separate  the  inhomogeneity  of  the 
refractive  index  from  the  extinction  coefficient.  Of  the  many 

computations  chat  have  been  performed  we  simply  show  two  of  the  most 
characteristic  ones. 

The  simplest  case  is  that  of  a  homogeneous  layer.  Taking  starting 
values  of  n  ■  2.3  and  k  -  10~',  the  optical  constants  were  calculated 
for  two  wavelengths  400  and  800  nm  and  the  indices  found  were  to  be 
within  0.03  2  of  the  initial  value.  As  foreseen,  the  efficiency  of  the 
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calculation  of  the  extinction  coefficient  was  not  as  good,  leading  to  a 
10  Z  relative  error. 

To  test  the  derivation  in  the  case  on  an  inhomogeneous  layer,  an 
arbitrary  curve  of  refractive  index  was  chosen  to  compute  the  evolution 
of  transmission  :  n(d)  »  2  +  0.3  exp  (-d/200)  where  d  is  the  thicltness  of 
the  layer.  The  extinction  coefficient  was  assumed  constant  equal  to 
10~*.  Table  1  shows  the  calculated  results  for  two  wavelengths.  The 
columns  An^/n^g  etc.  list  the  error  in  the  determination  expressed  as  a 


percentage 


TABLE  1 


Wavelength 


(nm) 

nin 

Anin/nin 

nout 

4nout/nout 

lc  i 

400 

2.297 

0.152 

2.012 

0.09  2 

1.03x10“* 

800 

2.283 

12 

2.025 

0.7  2 

0-1.7x10“* 

Results  of  simulation  calculations 


These  results  indicate  that  we  can  expect  acceptable  accuracy  in  the 
refractive  index  profile  using  this  technique.  It  is  important  to  note 
that  the  higher  accuracy  corresponds  to  the  shorter  wavelength.  This  is 
due  to  the  greater  number  of  extrema  at  the  shorter  wavelength,  which 
provides  more  Information  so  the  envelopes  calculated  have  a  more  exact 
position.  This  inaccuracy  of  the  envelopes  leads  to  questionable  values 
of  extinction  coefficient  for  the  longer  wavelength.  These  results  are 
rigorously  true  only  for  the  cases  considered  but  they  strongly  suggest 

what  we  can  generally  expect  from  this  method.  The  relative  accuracy 
of  the  refractive  index  depends  also  upon  the  quality  of  the 
transmittance  measurements:  it  ranges  from  AT/T  to  5AT/T  where  AT/T  is 
the  relative  accuracy  achieved  in  the  transmittance  measurements. 

IV.  Experimental  results 

The  optical  constant  determination  has  been  carried  out  for  titanium 
dioxide  layers  only.  We  will  present  the  results  obtained  with  two 
different  layers,  the  first  layer  fitting  the  model  used  In  this 
derivation  and  the  second  being  unstable. 

The  starting  material  used  in  our  experiments  was  Tit0,  evaporated  by 
electron  bombardment  onto  a  glass  substrate. 

For  the  first  layer  we  consider,  the  oxygen  partial  pressure  was 
4. Ax  10“'  mbar  and  the  chamber  temperature  varied  from  204°  to  227°C 
during  the  deposition. 


The  profile  of  refractive  index  and  of  extinction  coefficient  are  plotted 
in  Fig.  3  for  the  wavelength  678  an.  The  thickness  of  the  layer  is  670 
am,  the  inner  index  is  2.135  and  the  outermost  index  1.794.  The  layer  is 
inhomogeneous,  as  expected  for  an  oxide  layer,  and  slightly  absorbant. 

The  correlation  between  index  and  extinction  coefficient  can  be 
interpreted  as  an  increase  of  the  degree  of  oxidation  with  the  thickness 
of  the  layer  and  as  a  decrease  of  the  packing  density  due  to  a  conical 
form  of  the  columnar  structure.  It  is  difficult  to  choose  which  has  the 
primary  effect. 

Fig.  4  gives  the  dispersion  of  the  innermost  and  outermost  refractive 
indices.  A  high  dispersion  for  the  outermost  index  suggests  a  more 
absorbant  outer  part  of  the  layer. 

The  layer  appears  to  be  stable  and  therefore  our  results  indicate  that 
it  is  an  inhomogeneous  layer  showing  a  decrease  in  its  packing  density 
and  a  decrease  in  its  oxidation  degree  in  the  direction  of  growth. 

We  now  examine  an  unstable  layer.  This  layer  was  deposited  in  a  chamber 
at  a  temperature  of  260°C  in  the  presence  of  oxygen  at  a  partial 
pressure  of  1.3  10-'  mbar. 

This  atmosphere  was  intentionally  deficient  in  oxygen  compared  with  the 
usual  conditions  for  the  deposition  of  titanium  dioxide  and  we  expected 
some  absorption  in  the  layer.  We  are  not  disappointed  as  shown  Fig.  5 


but  furthermore  another  phenomenon  seems  to  occur.  The  curve  of  the 
extinction  coefficient  falls  with  thickness  towards  a  value  which  could 
even  become  negative.  This  is  almost  certainly  the  effect  of  instability 
in  the  oxidation  of  the  layer:  starting  with  a  high  deficiency  of  oxygen 
the  inner  part  of  the  film  is  gradually  oxidized  as  the  layer  grows. 
This  defeats  the  technique  presented  here  since  the  extrema  used  to 
calculate  the  envelopes  are  changing  with  the  oxidation  of  the  layer  and 
also  because  we  have  assumed  the  innermost  refractive  index  to  be 
constant.  The  real  curve  is  impossible  to  obtain  but  we  can  expect  that 
it  would  present  a  less  inhomogeneous  profile. 

Other  experiments  have  been  carried  out  to  understand  the  limitations  of 
this  method.  The  instability  of  the  layer  can  manifest  itself  by  the 
occurance  of  negative  extinction  coefficients,  by  thicknesses  varying 
with  wavelength  or  by  very  misshapen  profiles  of  index.  Nevertheless, 
the  application  of  this  technique  to  stable  layers  gives  very  interesting 
results.  These  can  give  information  on  the  structure  of  the  layer  in 
terms  of  packing  density  as  well  as  in  terms  of  degree  of  oxidation.  It 
is  also  important  to  know  that  layerB  can  be  unstable  and  to  be  able  to 
recognise  such  an  instability.  Further  work  is  required  in  this  area. 

We  hope  that  eventually  it  may  be  possible  to  distinguish  an  instability 
caused  by  a  structural  rearrangement  from  a  reoxidation  of  some 
Innermost  portions  of  the  film. 


V.  Conclusion 


This  mechod  has  been  developed  for  layers  presenting  a  high  index,  a 
small  homogeneity,  and  a  small  extinction  coefficient.  It  permits  the 
determination  of  the  profiles  of  the  optical  constants  and  the  dispersion 
of  index  in  vacuo  provided  the  assumption  of  stability  is  fulfilled.  It 
identifies  layers  that  are  unstable.  It  makes  possible  the  study  of  the 
variations  of  the  optical  constants  with  changes  in  the  deposition 
parameters. 
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FIGURE  CAPTIONS 


Plot  of  noisy  signal.  The  extrema  are  difficult  to  determine 
with  accuracy. 

After  filtering,  the  extrema  have  been  extracted  from  the  noise 
without  distorting  or  attenuating  the  transmission  curve. 

Profile  of  refractive  index  and  extinction  coefficient  for  a 
stable  titanium  dioxide  layer.  (Upper  curve  represents  n, 
lower  curve  k) . 

Dispersion  of  innermost  and  outermost  refractive  index  for  a 
stable  layer  of  titanium  dioxide. 

Example  of  result  given  by  the  method  when  applied  to  an  unstable 
layer.  Titanium  dioxide  layer  deposited  in  an  oxygen  deficient 
atmosphere. 


100H  FILTERED  TRANSMISSION 


WAVELENGTH  678  nm 
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of  resuit  given  by  Che  method  when  applied  to  an  unatable 
Titanium  dioxide  layer  deposited  in  an  oxyee„  deficient  ntmosphe 
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